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ABSTRACT 

We investigated the copper abundances for 64 late-type stars in the Galactic disk and halo with 
effective temperatures from 5400 K to 6700 K and [Fe/H] from —1.88 to —0.17. For the first time, 
the copper abundances are derived using both local thermodynamic equilibrium (LTE) and non-local 
thermodynamic equilibrium (non-LTE) calculations. High resolution (i? > 40,000), high signal-to- 
noise ratio {S/N > 100) spectra from the FOCES spectrograph are used. The atmospheric models are 
calculated based on the MAFAGS opacity sampling code. All the abundances are derived using the 
spectrum synthesis methods. Our results indicate that the non-LTE effects of copper are important 
for metal-poor stars, showing a departure of ^ 0.17 dex at the metallicity ^ —1.5. We also find 
that the copper abundances derived from non-LTE calculations are enhanced compared with those 
from LTE. The enhancements show clear dependence on the metallicity, which gradually increase with 
decreasing [Fe/H] for our program stars, leading to a flatter distribution of [Cu/Fe] with [Fe/H] than 
previous work. There is a hint that the thick- and thin-disk stars have different behaviors in [Gu/Fe], 
and a bending for disk stars may exist. 

Subject headings: Galaxy: evolution - line: formation -- line: profiles — stars: abundances - stars: 
late-type 


1. INTRODUCTION 

The investigation of the elemental abundance traces 
the evolution of the Galaxy. The chemical history of 
the Galaxy is dominated by the processes of nucleosyn¬ 
thesis in each generation of stars. Based on the ob¬ 
servational behaviors of elemental abundances in stars 
with different metallicities, one can not only look back 
the Galactic chemical enrichment history, but also con¬ 
strain the theoretical evolutionary models of our Galaxy. 
Copper is of particular interest among the iron-peak el¬ 
ements, because firstly it has unique evolutionary trend 
as the [Fe/H] varies from extremely metal-poor to the 
solar abundance, and secondly this element is thought 
to be synthesized by several possible nucleosynthesis sce¬ 
narios, yet the contributions of these scenarios are still 
in dispute. 

From the observational point of view, although it was 
iCohenl (|1978l . I1979LI1980D who gave the early glance at 
the Galactic copper abundances (they investigated 27 
red giants in 7 globular clusters), and who first sug¬ 
gested a decreasing trend of [Cu/Fe] with decreasing 
[Fe/H], the evolutiona ry trend of [Cu/Fe] was not es¬ 
tablished fir mly until iSneden et al~ (119911). In a se¬ 
ries of work (iGratton fc Snedenlll988l : lSneden fc Crocked 
119881: ISneden et al.l 119911) . the authors derived the cop¬ 
per abundances for a large sample of stars in the 
Galactic disk and halo. Their results clearly showed 
a linear trend of [Cu/Fe], which increases towards the 
higher metallicity. Th is tre nd was partly confirme d by 
iMishenina et al.l (|2002D and lSimmerer et (|2003D . the 
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former authors investigated an expanded sample of disk 
and halo stars with wide range metallicity, while the lat¬ 
ter ones managed to measure the copper abundances 
for 117 giants in 10 globular clusters. Furthermore, 
both work indicated a flat plateau of [Gu/Fe] at the 
metal-poor end, which is roughly [Cu/Fe] —0.75 at 
[Fe/H] < —1.5. Some contemporaneous studies on sin¬ 
gle ultra-metal-poor s tars seemed consistent with this 
flat distribution (^e.g.. iWestin et ahi 120001 : iCowan et al.l 
120021: ISneden et al.l 12(1(131). How ever, using the near- 
UV lines. iBihain et al.l (|2004D and lLai et al.l (|2008D sug- 
gested the pl a teau should be around —1 instead of —0.75. 
iReddv et al.l (|2003L 120061) investigated a large sample of 
disk stars, and they found little variation of [Cu/Fe] with 
[Fe/H] in the metallicity range [Fe/H] > —0.8, where the 
values of [Cu/Fe] show no evident divergence with re¬ 
spect to the solar copper abundance. Thus, a slight S- 
shape of [Cu/Fe] as a function of [Fe/H] c an be seen by 
overl apping the data from different authors ([Bihain et al.l 
12004 Figure 1). Although the stars investigated by a 
great quant ity of work followed the Galactic general S- 
shape (e.g.. iProc haska 6/^1120001: iShetron e et' al.ll200lL 
1200,81 : iCohen et al.l 120081 : IMishenina et al.l 120111) . several 
peculiar structures were still detected, such as U rsa Ma¬ 
jor moving Group (U MaG) (iCastro et alJll99^. Omega 
Cent auri (w Gen) ([Smith et al.l 120001: iCunha et al.l 
I2002D. Sagittarius dwarf sphe r oidal g alaxy (Sgr dSph) 

( McWilliam fc Smecker-Hanel 1200,51 : iMcWillian^^e^M 
20^), and the halo sub-population (iNissen fc Schuster 
20111 ). Their discordant trends of [Cu/Fe] may imply 
different chemical evolutionary histories with respect to 
our Milky Way. 

Theoretically, copper is thought to be produced in 
multiple astrophysical sites. The first one is the weak 
s—process which takes place in massive st ars during the 
heliu m- and carbon-shell burning stage (jSneden et al.l 
119911) . It is a secondary process that needs iron seeds 
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from previous generations of stars, result ing in the lin¬ 
ear d ependence of [Cu/Fe] on [Fe/H]. iBisterzo et alJ 
(|2004ll proposed a revised version of this scenario, sug¬ 
gesting sr—process dominated copper synthesis in mas¬ 
sive stars instead of classical s—process. Additionally, 
simulations of Galact i c che mical evolution (GCE) by 
iRomano fc Matteuccil (|2007ll and iRomano et alJ (120101 1 
also support t hat massive star s cont ribute most of cop¬ 
per. However. iMatteucci et al.l (1199311 showed the second 
possibility by fitting GCE models to interpret observa¬ 
tional data. They suggested that the main source of cop¬ 
per was type la supernovae (SNe la) instead of s—process 
(but they needed to increase the yields from SNe la by 
about an order of magnitude). This co nclusion was sub¬ 
sequen tly s upported by the work from iMishenina et all 
(j2002ll and iSimmerer et al.l (120031 1. The third mecha¬ 
nism comes from the constraint of the [Cu/Fe] plateau 
at low metallicity, which requires the contributions 
from primary explosive n ucleosynthesis in ty pe II su¬ 
pernovae (SNe II) (e.g., iTimmes et al.l 1199^ . How¬ 
ever, the pr oductions calculated from S Ne H are model- 
dependent (iWoo slev fc Weaver! 119951: iKobavashi et al.l 


120001 : iNomoto et al.1120(101 : IRomano et al.[ 20101 1. result- 
ing in discordant predicted fractions. Besides, large 
scatters in observational results also play a relevan t 
role (see the examples given bv iPignatari et all 1201011 . 
The last source was the main s—process operating in 
the low- and intermediate-mass AGB stars, which was 
thought to contribute only about 5% of the solar copper 
(|Travaglio et al.l [200411 . 

So far, a number of analyses on copper abundance have 
been presented, covering the whole range of metallicity, 
from solar to the most extreme metal-poor stars, but 
none of them is carried out with the non-local thermody¬ 
namic equilibrium (non-LTE) calculations. This is partly 
because non-LTE calculations need reliable atomic model 
and detailed statistical equilibrium calculations, whereas 
neither of them is a simple job for copper. Despite the 
difficulties, there are at least two principal reasons that 
we have to put non-LTE analysis into perspective. First, 
recent studies on non-LTE effects have demonstrated 
that the non-LTE correcti on is large for some elements 
in metal-poor stars (e.g., iBaumuel l er fc Ge h renI [r997l: 
Baumueller et al.lll994lZhao et aTIIT998l:lZhao fc GehrenI 


20001: Gehren et al. 12004 iBergeinann fc GehrenI 120081 
Shi et al.ll2004l20'(i^ . Moreover, only a few lines can be 


used to perform copper abundance analysis when dealing 
with the metal-poor stars with relatively higher temper¬ 
ature, as there is little neutral copper in the atmospheres 
of such stars. The analysis has to rely on the two res¬ 
onance lines at 3247 A and 3273 A, both of which ma y 
suffer large n on-LTE effe cts (IR oederer et al.ll2 012L 120141) . 
iBihain et al.l (|2004f) and iBonifacio et al.l ( 201011 also re¬ 
ported that the abundances derived from these two lines 
are inconsistent with those derived from other optical 
Gu I lines. 

In this paper, we aim at exploring the copper abun¬ 
dances for the sample stars in the metallicity range 
— 1.88 < [Fe/H] < —0.17 with a complete spectrum 
synthesis method based on level population calculated 
from the statistical equilibrium equations. In Section 
2, we will briefly introduce some key information about 
the observations. The fundamental work for abundance 


analysis, i.e., atmospheric model, stellar parameters, and 
atomic line data are presented in Section 3. The 4th 
section will provide details of non-LTE calculations in¬ 
cluding the atomic model. In Section 5 we will show 
the final results and error analysis. The discussions and 
conclusions are presented in Section 6 and Section 7, re¬ 
spectively. 

2. OBSERVATIONS 

The sample stars investi gated in this p aper h ave al - 
ready been discussed by IGehren et al.l (|2004 1200611 . 
Here, we simply list the key features of the sample and 
the observations. More details can be found in aforemen¬ 
tioned papers. 

- The observations were carried out on the 2.2 m tele¬ 
scope located at Calar Alto Observatory from the 
year of 1999 to 2003. The FOCES echelle spec¬ 
trograph was used to obtain the high-resolution 
spectra, providing 97 spectral orders in total that 
started at 3700 A and ended at 9800 A. 

“ The detector was a GGD chip with 2048 x 2048 
pixels, and the size for each pixel is 24 /im. A two- 
pixel bin results in a ^ 40, 000 resolution power 
(R). 

- The total exposure time is divided into more than 
two exposures. The final combined spectra show a 
signal-to-noise ratio (S/N) that is higher than 100. 

- The spectra were reduced from the program de¬ 
signe d for FOCES spectrograph (I Pfeiffer et al.l 
1199811 . which worked under the IDL environment. 
Cosmic rays and bad pixels were removed by care¬ 
ful comparisons of the exposures from the same ob¬ 
ject. The instrumental response and background 
scatter light were also considered during the data 
reducing. 

3. FUNDAMENTAL WORK FOR ABUNDANCE ANALYSIS 

In this section, we briefly describe the methods and 
key information our studies based on, including the at¬ 
mospheric model, the determinations of stellar parame¬ 
ters, the atomic line data, and the kinematic properties 
of our program stars. 

3.1. Model Atmosphere 

Stellar atmospheric model is the foundation of the 
spectrum synthesis. Our work adopted the MAFAGS 
op acity sa r npling (OS) model. This code was developed 
by iGrupd (|2004 and updated by iGrupp et al.l (1200911 
with the new iron atomic data computed bv iKurucz I 
(|2009fl . MAFAGS OS code describes a one-dimensional 
plane-parallel model with 80 layers overall in the hy¬ 
drostatic equilibrium state. The chemical homogeneity 
and local thermal equilibrium is assumed throughout the 
atmosphere. This m odel atmosphere was also applied 
in ou r previous work (iMashonkina et aLll2011l : iShi et ahl 

i2?fn . 

The comparison between MAFAGS and M ARGS OS 
model has been performed bv iShi et al.l (120141 1. Although 
differences of the temperatures between those two models 
exist outside logr ~ —4 and inside logr cs 0.3, both of 
these regions barely influence the synthetic copper spec¬ 
tral line profile. 































































































Non-LTE analysis of neutral copper 


3 


TABLE 1 

Atomic data of copper lines used in this work 


•^air 

(A) 

Transition 

-^low 

(eV) 

logs/ 

logCe 

5105.541 

4 -" "- 05/2 

- 4p"P3°/2 

1.389 

-1.64 

-31.67 

5218.202 

4p"R3/2 

- 4<i"D5/2 

3.817 

+0.28 

-30.57 

5782.132 


-4p=G% 

1.642 

-1.89 

-31.66 


Notes. The log 5 / values were rectified from the non-LTE solar 
spectrum fitting, and the v an der Waals damping constan ts (logCe) 
were calculated according to lAnstee Sz O’Maral il991lfl995l) 

3.2. Stellar Parameters 

For all of our program sta rs, we directly adop t ed th e 
stellar parameters derived bv iGehren et alJ (j2004lW)^ . 
By fitting the theoretical profiles to the observational 
data, two Balmer lines were used to derive the effective 
temperatures (Teff)- The hydrogen broadening theory 
involved in the the oretical profile calculation was from 
lAli fc GriemI (1196611 . The surface gravities (log 5) were 
obtained by [g] = [Mass/Luminosity] + 4[reff], where 
the stellar mass and lumin osity were evaluated with the 
help of evolutionary tracks (|VandenBerg et ani2000ll and 
HIPPARCOS parallaxes. In addition, the microturbu¬ 
lence velocities (^) were determined simultaneously with 
the metallicities ([Fe/H]): the [Fe/H] determined from 
Fe II line is supposed to be independent of the equiva¬ 
lent width. The final uncertainties in Tea, logg, [Fe/H], 
and ^ were estimated to be ±80 K, ±0.05 dex, ±0.05 
dex, and ±0.1 km s“^, respectively. 

3.3. Atomic Line Data 

We used a set of calibrated log gf values, each of which 
could reproduce the so lar coppe r abun dance indepen¬ 
dently, as presented bv IShi et al.l (|2014f) . Furthermore, 
the van der Waals damping c onstants (logCfi) for Cu I 
were calculated according to lAnstee fc O’Mar^ (I1991L 
Il995f) . Five Cu I lines can be seen in our FOCES spectra. 
For the evaluation of copper abundance, 5220.070 A and 
5700.240 A are not good indicators, for both of them are 
badly blended in most of our program stars. As a result, 
the rest three lines are applied to our research, which are 
5105.541 A, 5218.202 A, and 5782.132 A. The atomic line 
data of the three lines are listed in Table [T] 

3.4. Population and Kinematic Properties 

Similar to the stellar p arameters, we also a d opted the 
population identified by IGehren et al.l (|2004 I2006f) for 
our program stars, which was based on the kinematic 
features, stellar ages, [Al/Mg] and [Mg/Fe] ratios. Most 
of the stars were classified as thin-disk, thick-disk, and 
halo population, whereas the rest were the stars with 
peculiarities. 

4. NON-LTE CALCULATIONS 

LTE assumption provides us a simple way to calculate 
the population of each energy level and number densities 
of different ionization stages for a given element, while 
non-LTE calculations require to solve the detailed sta¬ 
tistical equilibrium equations. Thus, a reliable atomic 
model of copper is indispensable. 



Ug T^hkki 

Fig. 1.— The departure coefficients {bi) for selected energy levels 
(listed in the figure) as a function of continuum optical depth at 
5000 A for the model atmosphere of HD 59984. The coll ision with 
neutra l hydrogen was scaled by a factor of 0.1 in light of IShi et al.l 
ll20ll'l . 

The atomic m odel of copper h as been described in the 
previous paper (|Shi et al.l 12014 and the Grotrian dia¬ 
gram of the model can also be seen there. We mod¬ 
eled the copper atom with 17 orbits, 97 energy levels (96 
states for Cu I and the ground state for Cu II) and 1089 
transitions, and the fine structure for the levels with low 
excitation energy was also included. The atomic data of 
such complex structure are obtained from both labora¬ 
tory measurements (NISlQ databas e. iSiigar fc Mnsgro^ 
1 19901 ) and theoretical calculations (|Liu et al.l 120141) . Ti 
addition, the excitation and ionization caused by inelas¬ 
tic collisions were also considered. The data for col- 
lisio ns with neutral hydrogen were obta ined based on 
the iDrawinI formula (iDrawinI I1968L I1969D presented by 
iSteenbock fc Holwegen i 19841) . and we decreased the col- 
lisional rates by an o rder of magnitu de (5f= 0.1) un¬ 
der the suggestion of iShi etlin (12014 . The excitation 
and ionization caused by inelastic collisions with elec¬ 
trons are calculated according to a numb e r of theoretica l 
work ([van Regemorteiiri962l : IAllen~lll973l: iSeaton II1962D 
We u sed a revised DETAIL program ([Butler fc Giddinra 
1 19851 ) with accelerated lambda iteration method to per¬ 
form the statistical equilibrium calculations. 

In Figure (U we present how the departure coefficients 
{bi = of the selected levels vary with 

the continuum optical depth at 5000 A (logrsooo) for the 
model atmosphere of HD 59984 , where bi is the ratio of 
the number density of non-LTE to that of 

LTE (n^""^^). HD 59984 is a typical star with moderate 
temperature and metallicity among our sample, which is 
randomly selected as an example for the convenience of 
discussion. The departure coefficients for Cu I important 
levels and Cu II ground state are shown in the figure. It 
also shows that the number densities of these levels begin 
to underpopulate outside the layers with logrsooo ~ 0.5 
due to the overionization. 

5. RESULTS 

5.1. Spectral Line Synthesis 

^ http://www.physics.nist.gov/ 
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Fig. 2.— The synthetic profile of the Cu I 5105 A line for 
HD 59984. The observed spectrum and theoretical synthesis are 
represented by filled circles and solid line, respectively. 

In our analysis, the transitions of the hy perfine struc - 
ture (HFS) were calculated according to iBiehll (Il976f l 
with RS coupling method. For the lines we used, if the 
intervals of HFS components are within 1 mA, we com¬ 
bined them together. The adopted solar copper abun¬ 
dance is the value derived from the me teorites, which is 
log£ 0 (Cu) = 4.25 (|Lodders et al.l[200^ . and the oscilla¬ 
tor strengths were also calibrated based on this value as 
described in aforementioned section. Additionally, the 
ratio between two copper isotopes (^^Cu and ^ ^Cu) was 
assumed to be 0.69 : 0.31 (|AsDlund et al.ll200^ . An IDL 
based program SIU was used to perform the line forma¬ 
tion in our abu n dance determinations, which was devel¬ 
oped bv IReetz! (1199111 . In our analysis, the broadenings 
caused by the macroturbulence, rotation, and instrument 
were treated as one single Gauss broadening factor, being 
convolved with the synthetic spectra to fit the observed 
line profile. The comparison between the synthetic and 
observed line profile at 5105 A for HD 59984 is shown in 
Figure [5] as an example. The observed spectrum and 
theoretical synthesis are represented by filled circles and 
solid line, respectively. The uncertainty of our line profile 
synthesis is less than 0.02 dex. 

5.2. Copper Abundances and Error Analysis 

The copper abundances are derived successfully for 
60 stars in our sample, and the selected Cu I lines are 
too weak to rely on for the rest four stars, which are 
HD 241253, HD 233511, G 119-32, and BD 4-20° 2594. 
The derived abundances with both LTE and non-LTE 
calculations are listed in Table [H where the results of 
each individual line are also presented. The final abun¬ 
dance of each star is given by calculating the arithmetical 
mean value of every line used in the analysis. 

The trend of [Cu/Ee] to [Fe/H] for our program stars 
is shown in Figure |51 The errors in the figure are evalu¬ 
ated by computing the standard deviations of the abun¬ 
dances derived from different spectral lines. Our results 
do not show large abundance discrepancy between differ¬ 
ent lines. We present the errors as a function of [Fe/H] 
in Figure |4| In LTE calculations, the standard devia¬ 
tions vary from 0.01 to 0 . 11 , slightly larger than those 


in non-LTE, which are between 0.01 and 0.08. Further¬ 
more, the LTE standard deviations become large at the 
metal-poor end, while the non-LTE ones remains stable. 
Both of them have a mean value that are around 0.04. 
The errors caused by the uncertainties of the stellar pa¬ 
rameters are estimated for HD 59984, and the resulting 
effects in [Cu/Ee] are ±0.07, < 0.01, ±0.05, and < 0.01 
dex for the typical uncertainties in logg, [Ee/H], and 

respectively. 

The differences in [Cu/Ee] between non-LTE and LTE 
for our program stars as a function of metallicity, effective 
temperature, and surface gravity are plotted in Eigure[SJ 
The non-LTE departure shows clear dependence on the 
metallicity (see Eigure[5^), which is gradually increases 
as [Ee/H] decreases in our program stars. 

5.3. Non-LTE Effects 

The final results in Table [2] show that the abundances 
derived from non-LTE calculations are larger than those 
from LTE ones for our program stars. The non-LTE cor¬ 
rection for individual spectral line can reach ^ ± 0.20 dex 
for the stars with [Fe/H] ^ —1.5. The increasing of the 
copper abundance in non-LTE calculation is the conse¬ 
quence of the underpopulation of the lower energy levels 
(4s^^D and 4p^P°), as shown in EigureHj This devia¬ 
tion will lead to a underestimation of copper abundance 
by LTE analysis. 

For each Cu I line, the non-LTE effect is different, re¬ 
flecting the properties of corresponding energy levels and 
associated transitions. The lines of 5105 A and 5782 A ex¬ 
hibit larger non-LTE effects compared with the weaker 
5218 A. Taking HD 59984 as an example again, the non- 
LTE corrections for the lines of 5105 A and 5782 A are 
0.10 and 0.09 dex, respectively, while it is 0.06 dex for 
the weaker 5218 A line. On one hand, despite the 5218 A 
line suffers less non-LTE effects than the other ones, it 
is usually blended by the line of 5217 A (Fe II line), thus 
it is not a satisfactory indicator for the stars with low 
or moderate metallicity. On the other hand, since most 
abundance analyses are based on the two copper strong 
lines, one need to be aware of the non-LTE departures, 
especially for the metal-poor stars. 

5.4. Comparison with Other Work 

Numbers of studies for copper abundance have been 
carried out so far, allowing us to compare our results 
with those derived by several different groups based on 
LTE calculations. Comparing with the work published 
in recent years, we found 47 stars are in common (11 
of them are studied at least twice in different papers). 
In Eigure | 6 l we compare our results derived from LTE 
calculations with those in the other work. In general, 
the LTE [Cu/Fe] values determined in our work are con¬ 
sistent with the different studies presented here, and no 
clear systematic deviation is found. We briefly discussed 
the details of the comparisons and the possible reasons 
o f the large scatter s. _ 

\Mishenina et al\ K200A \20lA }: iMishenina et al.l stud¬ 
ied cop per a bundances for 90 metal-poor stars in the 
year of 120021 and later they investigated 172 F to K 
dwarf stars and derived their copper abundances in l 201 lL 
All their observations were performed with R = 42,000 
and S/N > 100. The copper abundances were derived 
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- 1.5 - 1.0 - 0.5 0.0 

[Fe / H] 


Fig. 3.— Abundance ratios [Cu/Fe] as a function of [Fe/H] for our program stars, where non-LTE and LTE results are represented by 
filled and open symbols, respectively. Furthermore, symbols with different shapes represent stars from different populations, which are: 
triangle - thin disk, circle - thick disk, star - halo, square - objects with peculiarities. The errors are evaluated by computing the standard 
deviations of the abundances derived from different spectral lines. 


0.0 -0.5 -1.0 -1.5 

[Fe/H] 

Fig. 4.— The distribution of standard deviation as a function 
of metallicity, where filled and open circles represent non-LTE and 
LTE calculations, respectively. The stars whose copper abundances 
were derived by only one single line were not plotted in the figure. 

from the 5105 A, 5218 A, and 5782 A lines, which were 
the same as those in our studies. The oscillator strengths 



adopt ed in their studies were from iGurtovenko fc Kostvfl 
(|1989r) . We have 13 stars in common, and four of them 
are investigated in both of their papers. Since the derived 
[Cu/Fe] of these four stars varied a little bit for their two 
pape rs, we therefore adopted the latest values (derived in 
1201 If) to perform the comparison. The resulting average 
difference of the 13 stars is A[Cu/Fe] = —0.01±0.12. The 
main contribution of the relatively large scatter is from 
four stars, namely HD 22879, HD 101177, HD 108076, 
and HD 218209, as the [Cu/Fe] values derived from our 
work and theirs differ by more than 0.1 dex for those 
stars. This is mainly due to the differences of the stellar 
p arameters ad o pted i n the two studies. 

\Reddv et al\ i2003l . \200A) : The authors performed an 
abundance analysis on a large sample of F and G dwarfs, 
containing thin-/thick-disk and halo stars. The copper 
abundances were determined from the 5105 A, 5218 A, 
and 5220 A lines. Their log 5 / values are similar to 
ours. From their samples, 21 stars are found in com¬ 
mon with ours, and the mean difference is A[Cu/Fe] = 
0.01 ± 0.10. Most of the stars show no large deviation 
except HD 200580 and HD 107582. The [Fe/H] values for 
these two stars adopted in their studies are different from 
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Fig. 6.— Comparison of derived [Cu/Fe] in LTE analy¬ 
sis for the stars in common with th e literature, includin g 
MKS02 l[Mishenin^eta2j2Q03) , MGBjd_J[Mjahenimrjt an i201l| i. 
RTL03 [200^, RLA06 l|Reddy_e^a|J[2b03) , ABL0 4 

(lAllendePrietoer^ I 2 SMI . and NSll (lNTSetr7rSchustCTl2011l i. 
Corresponding symbols are annotated in the figure. 

derive copper abundances. The log gf values adopted 
in their analyses are similar to those used in ours. Our 
results are quite consistent with theirs. The average dif¬ 
ference is 0.00 ± 0.05 for the 11 stars in common. 

6. DISCUSSION 



Fig. 5.— The differences in [Cu/Fe] between non-LTE and LTE 
for our program stars as a function of metallicity (a), effective 
temperature (b), and surface gravity (c). 

those in ours. Taking [Fe/H] variances into account, the 
d eviation can be perfec t ly rem oved. 

\Allende Prieto etaU H200A) : These authors analyzed 
a complete and comprehensive sample of 118 stars with 
absolute magnitude brighter than 6.5 and distance less 
than 14.5 pc from the Sun. The spectral were obtained 
with R = 50, 000 and S/N > 150. Compared with their 
sample, we have only one star in common. The difference 
o f the [Cu/Fe] value i s —0.0 7. _ 

\Nissen & Schust^ 11201 A) : iNissen fc Schusteil carried 
out a series of elemental abundance analyses on a sample 
of 94 dwarfs, most of which were identified as halo stars. 
They found the halo stars in the solar neighborhood fall 
into two distinct populations that can be separated by 
[a/Fe]. They used 5105 A, 5218 A, and 5782 A lines to 


6.1. The Evolutionary Trend of [Cu/Fe] and 
Nucleosynthesis in the Galaxy 

The observational trend of [X/Fe] as a function of 
[Fe/H] is a powerful tool to reveal the origins of elements 
and constrain the Galactic chemical evolutional model. 
The trend of copper for our Galaxy has been investigated 
in many papers, but all of the calculations are under LTE 
assumptions. 

Figure [3] displays the results of [Cu/Fe] in our program 
stars as a function of [Fe/H] for both LTE and non-LTE 
calculations. The [Cu/Fe] trend is similar to the ear¬ 
lier work for the LTE results, while it is not the case 
for non-LTE ones. In order to show the features more 
clearly, we present the average [Cu/Fe] values for each 
0.1 dex bin of [Fe/H] (A[Fe/H] = 0.1 bin) in Figure [7] 
for non-LTE results. Stars from different populations 
are averaged separately, and the peculiar stars are not 
included. The error bars represent the abundance dis¬ 
persions of the corresponding bins (no bar is plotted if 
there is only one or two stars in the bin). A flat distri¬ 
bution of [Cu/Fe] can be seen for our non-LTE results 
in the range of —1.5 < [Fe/H] < —1.0, which is different 
from that revealed by the previous work (a linear in¬ 
crease from [Fe/H] = —1.5 to —1.0). However, it should 
be noted that there are only four stars in that region, 
and more data are needed to confirm this trend. Even 
though there may exist bias, the [Cu/Fe] trend derived 
from non-LTE is still much flatter than that from LTE at 
— 1.5 < [Fe/H] < —1.0. The [Cu/Fe] trend for d isk stars 
is not discussed often. iProchaska et al.l (|2000f l studied 
the copper abundances for 10 thick-disk stars, and a su¬ 
persolar [Cu/Fe] can be see n in their results at [Fe/H 
~ —0.4. iReddv et al.l (1200611 suggested that the [Cu/Fe 
seemed slightly greater for thick-disk stars than that for 
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Fig. 7.— The trend of average [Cu/Fe] for 0.1 dex bin of [Fe/H], 
only non-LTE results are presented. Stars from different popula¬ 
tions are averaged separately, and the peculiar stars are not in¬ 
cluded. The error bars represent the abundance dispersions of the 
corresponding bins. 


thin-disk stars. Our results indicate that the [Cu/Fe] 
gradually increases with the increasing [Fe/H] for the 
thick-disk stars, while most of the thin-disk stars have a 
solar [Cu/Fe] value in the overlapping region. Thus, the 
thick-disk population appears a slight overabundance of 
[Cu/Fe] at -0.7 < [Fe/H] < -0.4. 


Several g roups 


Timmes et al 


le-g-, 


IMatteucci et al.l 


]Mishenina^£^al 


Kobavashi et a.1.1 120001 iRomano fc Matteuccil 


1993; 


2002 ; 


2007; 


Romano et al.l 2010l) have modelled the Galactic chem¬ 


ical evolution of copper. The most difficult part in the 
modeling is to approach satisfactory copper abundances 
at both metal-poor end and solar- metallicity: normal 
type H supernovae yields from iWooslev fc Weav^ 
case B give a good approaching to the ob¬ 
served [Cu/Fe] trend for [Fe/H] > —2, but lead to 
overabundant f or copper at lower m etallicity, while 
the yields from iKobavashi et al.l (1200611 SNe H mixed 
with hypernovae give a better htting in the metal-poor 
end but fail to repr oduce the observatio nal trend at 
the solar-metallicity (|Romano et al.l 1201011 . The ways 
to solve this inconsistency are so limited because: (1) 
the co ntribution from s—pro cess cannot be changed 
freely (|Matteucci et al.l 1199311 and (2) there are still 
large uncertainties in modeling the yields from low- and 
intermedi ate-mass AGB stars or supe rn ovae explosion 
event s ([Romano fc Matteuccil 120071 : IRomano et al.l 
[Moll . Since the departures from LTE for Cu I show a 
clear dependence on metallicity, the copper abundances 
of the very/extremely metal-poor stars are expected to 
increase towards lower metallicity ([Roederer et al.l 12014 
Shi et al. 2014, in preparation). Consequently, before 
modifying GCE models, a firm and reliable observational 
trend of [Cu/Fe] should be established first. 


feature can be seen in F igure [7] as well . The term ‘b end¬ 
ing’ was introduced b y iBisterzo et al.l (12004 1200511 and 
IRomano fc Matteuccil ([200711 . Those authors superim¬ 
posed the [Cu/Fe] data from different work, and a bend¬ 
ing for disk stars appeared in the overlapping region, 
however, the possibility that the bending may be caused 
by the systematic offset between different studies could 
not be ruled out. The similar feature can be found for 
our non-LTE results, and it should be noted that we do 
not have enough data at [Fe/H] > —0.3, the decline in 
this region is mainly produced by three stars, and as a 
result, more data are needed to draw a hrm conclusion. 

One of the explanations for the bending might be 
SNe la. However, the typical timescale for the chemi¬ 
cal enrichment f rom SNe la in the Milky W ay is model- 
dependent le.g.. IMatteucci fc Recchiir2nnill and thus not 
strictly constrained. 

7. CONCLUSIONS 

We have investigated the copper abundances for 64 
late-type stars with effective temperatures from 5400 
K to 6700 K and [Fe/H] between —1.88 and —0.17. 
The non-LTE statistical equilibrium calculations are per¬ 
formed to derive the copper abundances. Based on our 
results, we come to the following conclusions. 

1. The non-LTE effects are signihcantly large for cop¬ 
per, and they differ from line to line. Two Cu I 
lines, 5105 A and 5782 A are more sensitive to non- 
LTE effects. The non-LTE corrections for these 
lines can reach ^ 0.2 dex at the metallicity [Fe/H] 
^ —1.5. The weaker line 5218 A shows relatively 
small but still not negligible non-LTE effects, with 
a maximum departure of 0.13 dex in our sample. 

2. The copper abundances are underestimated in LTE 
calculations. Taking non-LTE effects into account, 
the copper abundances increase for all of our pro¬ 
gram stars. 

3. The non-LTE effects clearly show dependence 
on metallicity, and they increase with decreasing 
[Fe/H]. 

4. Our non-LTE results show that there might be a 
Cu/Fe] plateau in the metallicity range —1.5 < 
Fe/H] < —1.0, however, it need more data to con¬ 
firm this result. There is a hint that the thick- and 
thin-disk stars have different behaviors in [Cu/Fe], 
and a bending for disk stars may exist. 

Further non-LTE studies to copper abundance in var¬ 
ious environments with wider metallicity range need to 
be carried out. They are essential to the understanding 
of the copper origination and evolution. 


6.2. the Bending 

The non-LTE results for disk stars suggest that there 
may be a bending-like feature at [Fe/H] > —1.0 (Fig¬ 
ure El), where the [Cu/Fe] trend goes up at first, but 
slightly decreases as the [Fe/H] becomes higher. The 
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TABLE 2 

Copper abundances of our program stars 


Star 

Te« 

log 9 

[Fe/H] 


5105A 

5218A 

5782A 

[Cu/Fe] : 

POP LinFor 

HD 17948 

6325 

4.13 

-0.35 

1.90 

-0.25 

-0.19 

-0.16 

-0.20 

± 

0.05 

D 

L 






-0.18 

-0.13 

-0.09 

-0.13 

± 

0.05 


N 

HD 22309 

5900 

4.29 

-0.31 

1.30 

0.03 

-0.04 

-0.03 

-0.01 

± 

0.04 

D 

L 






0.07 

-0.02 

0.01 

0.02 

± 

0.05 


N 

HD 22879 

5775 

4.26 

-0.83 

1.10 

-0.25 

-0.13 

-0.15 

-0.18 

± 

0.06 

T 

L 






-0.17 

-0.10 

-0.08 

-0.12 

± 

0.05 


N 

HD 30649 

5765 

4.26 

-0.58 

1.10 

0.03 

0.03 

0.02 

0.03 

± 

0.01 

T 

L 






0.09 

0.05 

0.07 

0.07 

± 

0.02 


N 

HD 243357 

5675 

4.38 

-0.59 

1.10 

-0.04 

0.05 

-0.02 

0.00 

± 

0.05 

T 

L 






0.01 

0.07 

0.03 

0.04 

± 

0.03 


N 

HD 36283 

5475 

4.28 

-0.41 

0.80 

0.01 

0.13 

0.05 

0.06 

± 

0.06 

T 

L 






0.04 

0.13 

0.08 

0.08 

± 

0.05 


N 

G 99 - 21 

5525 

4.30 

-0.63 

1.00 

-0.13 

0.01 

-0.10 

-0.07 

± 

0.07 

T 

L 






-0.09 

0.03 

-0.06 

-0.04 

± 

0.06 


N 

HD 250792 

5600 

4.32 

-1.02 

1.10 

-0.38 

-0.40 

-0.38 

-0.39 

± 

0.01 

H 

L 






-0.30 

-0.38 

-0.31 

-0.33 

± 

0.04 


N 

HD 46341 

5880 

4.36 

-0.58 

1.80 

-0.05 

-0.07 

-0.08 

-0.07 

± 

0.02 

T 

L 






0.01 

-0.04 

-0.02 

-0.02 

± 

0.03 


N 

HD 56513 

5630 

4.53 

-0.45 

1.20 

-0.06 

-0.01 

0.01 

-0.02 

± 

0.04 

D 

L 






-0.03 

0.00 

0.04 

0.00 

± 

0.04 


N 

HD 58551 

6190 

4.23 

-0.53 

1.80 

-0.13 

-0.07 

-0.02 

-0.07 

± 

0.06 

D 

L 






-0.05 

-0.01 

0.06 

0.00 

± 

0.06 


N 

HD 59374 

5840 

4.37 

-0.83 

1.40 

-0.16 

-0.17 

-0.08 

-0.14 

± 

0.05 

T 

L 






-0.08 

-0.13 

0.00 

-0.07 

± 

0.07 


N 

HD 59984 

5925 

3.94 

-0.74 

1.20 

-0.12 

-0.06 

-0.08 

-0.09 

± 

0.03 

T 

L 






-0.02 

0.00 

0.01 

0.00 

± 

0.02 


N 

HD 60319 

5875 

4.17 

-0.82 

1.40 

-0.15 

-0.10 

-0.15 

-0.13 

± 

0.03 

T 

L 






-0.06 

-0.05 

-0.06 

-0.06 

± 

0.01 


N 

G 235 - 45 

5500 

4.25 

-0.59 

1.10 

0.01 

-0.01 

0.05 

0.02 

± 

0.03 

T 

L 






0.05 

0.00 

0.09 

0.05 

± 

0.05 


N 

HD 88446 

5915 

4.03 

-0.44 

1.60 

-0.01 

-0.08 

-0.14 

-0.08 

± 

0.07 

D 

L 






-0.04 

-0.04 

-0.07 

-0.05 

± 

0.02 


N 

HD 88725 

5665 

4.35 

-0.70 

1.20 

-0.01 

0.03 

0.03 

0.02 

± 

0.02 

T 

L 






0.05 

0.06 

0.09 

0.07 

± 

0.02 


N 

HD 91784 

5890 

4.47 

-0.33 

1.30 

-0.01 

0.06 

0.06 

0.04 

± 

0.04 

D 

L 






0.03 

0.08 

0.10 

0.07 

± 

0.04 


N 

HD 94028 

5925 

4.19 

-1.51 

1.50 

-0.57 

-0.41 


-0.49 

± 

0.11 

H 

L 






-0.37 

-0.28 


-0.32 

± 

0.06 


N 

HD 96094 

5900 

4.01 

-0.46 

1.70 

-0.05 

-0.06 

0.01 

-0.03 

± 

0.04 

D 

L 






0.02 

-0.02 

0.07 

0.02 

± 

0.05 


N 

HD 97855A 

6240 

4.13 

-0.44 

1.80 

-0.22 

-0.15 

-0.11 

-0.16 

± 

0.06 

D 

L 






-0.14 

-0.09 

-0.03 

-0.09 

± 

0.06 


N 

HD 101177 

5890 

4.30 

-0.47 

1.80 

0.09 

0.04 

0.14 

0.09 

± 

0.05 

D 

L 






0.16 

0.07 

0.20 

0.14 

± 

0.07 


N 

HD 104056 

5875 

4.31 

-0.41 

1.30 

-0.02 

0.00 

0.02 

0.00 

± 

0.02 

D 

L 






0.03 

0.02 

0.06 

0.04 

± 

0.02 


N 

HD 107582 

5565 

4.34 

-0.61 

1.00 

-0.10 

-0.09 

-0.15 

-0.11 

± 

0.03 

T 

L 






-0.06 

-0.07 

-0.11 

-0.08 

± 

0.03 


N 

HD 108076 

5725 

4.44 

-0.73 

1.20 

-0.20 

-0.17 

-0.21 

-0.19 

± 

0.02 

D 

L 






-0.14 

-0.15 

-0.16 

-0.15 

± 

0.01 


N 

HD 114606 

5610 

4.28 

-0.57 

1.20 

0.04 

0.06 

0.02 

0.04 

± 

0.02 

T 

L 






0.09 

0.08 

0.06 

0.08 

± 

0.02 


N 

HD 118659 

5510 

4.36 

-0.60 

1.00 

-0.03 

0.01 

0.02 

0.00 

± 

0.03 

T 

L 






0.02 

0.02 

0.06 

0.03 

± 

0.02 


N 

HD 119288 

6420 

4.13 

-0.17 

1.90 

-0.34 


-0.26 

-0.30 

± 

0.06 

D 

L 






-0.28 


-0.20 

-0.24 

± 

0.06 


N 


Notes. Both LTE and non-LTE copper abundances (for each star, first and second row, respectively) of ou r program stars are listed in 
column 9. The abundances given here are the relative values to [Fe/H] derived from Fe II lines (Section l3.2ll . Column 6, 7, and 8 are the 
abundances derived from corresponding Cu I lines. The stellar parameters and population assignments are also shown in the table. The 
characters ‘D’, ‘T’, ‘H’ and in the ’POP’ column represent thin-disk, thick-disk, halo, and peculiar stars, respectively. The rightmost 
column indicates the line formation scenario for each star, where ’L’ represents LTE line formation and ’N’ represents non-LTE line 
formation. 
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TABLE 2 - continued 


Star 

T’eff 

logs 

[Fe/H] 


5105A 

5218A 

5782A 

[Cu/Fe] 

Pop 

LinFor 

HD 123710 

5790 

4.41 

-0.54 

1.40 

0.07 

0.08 

0.06 

0.07 + 0.01 

D 

L 






0.13 

0.11 

0.11 

0.12 + 0.01 


N 

HD 126512 

5825 

4.02 

-0.64 

1.60 

-0.03 

0.02 

0.06 

0.02 + 0.05 

T 

L 






0.05 

0.07 

0.14 

0.09 + 0.05 


N 

HD 134169 

5930 

3.98 

-0.86 

1.80 

0.01 

-0.04 


-0.02 + 0.04 

T 

L 






0.13 

0.04 


0.09 + 0.06 


N 

HD 142267 

5807 

4.42 

-0.46 

1.00 

0.01 

-0.06 

-0.02 

-0.02 + 0.04 

D 

L 






0.06 

-0.04 

0.02 

0.01 + 0.05 


N 

HD 144061 

5815 

4.44 

-0.31 

1.20 

0.00 

-0.06 

0.04 

-0.01 + 0.05 

D 

L 






0.04 

-0.05 

0.07 

0.02 + 0.06 


N 

HD 148816 

5880 

4.07 

-0.78 

1.20 

-0.04 

-0.07 

-0.02 

-0.04 + 0.03 

7 

L 






0.06 

-0.02 

0.07 

0.04 + 0.05 


N 

HD 149996 

5665 

4.09 

-0.52 

1.20 

0.00 

0.08 

0.00 

0.03 + 0.05 

T 

L 






0.06 

0.11 

0.05 

0.07 + 0.03 


N 

BD +68° 901 

5715 

4.51 

-0.25 

1.40 

-0.12 

-0.14 

-0.04 

-0.10 + 0.05 

D 

L 






-0.09 

-0.13 

-0.01 

-0.08 + 0.06 


N 

HD 157089 

5800 

4.06 

-0.59 

1.20 

-0.08 

-0.09 

-0.05 

-0.07 + 0.02 

T 

L 






-0.01 

-0.06 

0.01 

-0.02 + 0.04 


N 

HD 157466 

5990 

4.38 

-0.44 

1.10 

-0.20 

-0.22 

-0.15 

-0.19 + 0.04 

D 

L 






-0.14 

-0.20 

-0.10 

-0.15 + 0.05 


N 

HD 158226 

5805 

4.12 

-0.56 

1.10 

0.04 

-0.03 

0.04 

0.02 + 0.04 

T 

L 






0.11 

0.00 

0.10 

0.07 + 0.06 


N 

G 170 - 56 

6030 

4.31 

-0.79 

1.30 

-0.42 

-0.29 


-0.35 + 0.09 

7 

L 






-0.33 

-0.24 


-0.28 + 0.06 


N 

HD 160933 

5765 

3.85 

-0.27 

1.20 

-0.16 

-0.19 

-0.16 

-0.17 + 0.02 

D 

L 






-0.12 

-0.17 

-0.12 

-0.14 + 0.03 


N 

HD 160693 

5850 

4.31 

-0.60 

1.20 

0.05 

0.01 

0.08 

0.05 + 0.04 

7 

L 






0.12 

0.04 

0.14 

0.10 + 0.05 


N 

HD 170357 

5665 

4.07 

-0.50 

1.20 

-0.03 

-0.04 

0.03 

-0.01 + 0.04 

T 

L 






0.02 

-0.02 

0.08 

0.03 + 0.05 


N 

HD 171620 

6115 

4.20 

-0.50 

1.40 

-0.02 

-0.06 

0.01 

-0.02 + 0.04 

D 

L 






0.06 

-0.01 

0.08 

0.04 + 0.05 


N 

G 142 - 2 

5675 

4.48 

-0.75 

1.10 

-0.03 

-0.02 

-0.04 

-0.03 + 0.01 

T 

L 






0.03 

0.00 

0.01 

0.01 + 0.02 


N 

HD 182807 

6100 

4.21 

-0.33 

1.40 

-0.04 

-0.06 

0.00 

-0.03 + 0.03 

D 

L 






0.02 

-0.03 

0.06 

0.02 + 0.05 


N 

HD 184448 

5765 

4.16 

-0.43 

1.20 

0.21 

0.17 

0.22 

0.20 + 0.03 

T 

L 






0.27 

0.19 

0.27 

0.24 + 0.05 


N 

HD 186379 

5865 

3.93 

-0.41 

1.20 

-0.06 

-0.06 

-0.06 

-0.06 + 0.00 

D 

L 






0.00 

-0.03 

0.00 

-0.01 + 0.02 


N 

HD 198300 

5890 

4.31 

-0.60 

1.20 

0.09 

0.02 

0.06 

0.06 + 0.04 

T 

L 






0.16 

0.06 

0.12 

0.11 + 0.05 


N 

HD 200580 

5940 

3.96 

-0.82 

1.40 

0.16 

0.18 

0.28 

0.21 + 0.06 

7 

L 






0.28 

0.26 

0.39 

0.31 + 0.07 


N 

G 188 - 22 

6040 

4.37 

-1.25 

1.50 

-0.54 

-0.46 


-0.50 + 0.06 

H 

L 






-0.39 

-0.36 


-0.38 + 0.02 


N 

HD 201889 

5710 

4.05 

-0.78 

1.10 

-0.11 

-0.17 

-0.13 

-0.14 + 0.03 

T 

L 






-0.03 

-0.13 

-0.06 

-0.07 + 0.05 


N 

HD 204155 

5815 

4.09 

-0.66 

1.20 

-0.02 

0.00 

0.02 

0.00 + 0.02 

T 

L 






0.06 

0.04 

0.09 

0.06 + 0.03 


N 

HD 208906 

6025 

4.37 

-0.76 

1.40 

-0.07 

-0.09 

-0.02 

-0.06 + 0.04 

D 

L 






0.02 

-0.04 

0.06 

0.01 + 0.05 


N 

G 242 - 4 

5815 

4.31 

-1.10 

1.20 

-0.52 

-0.44 

-0.40 

-0.45 + 0.06 

H 

L 






-0.42 

-0.39 

-0.30 

-0.37 + 0.06 


N 

HD 215257 

6030 

4.28 

-0.58 

1.40 

-0.18 

-0.22 

-0.14 

-0.18 + 0.04 

D 

L 






-0.11 

-0.18 

-0.08 

-0.12 + 0.05 


N 

HD 218209 

5665 

4.40 

-0.60 

1.10 

0.14 

0.09 

0.17 

0.13 + 0.04 

T 

L 






0.20 

0.11 

0.22 

0.18 + 0.06 


N 

HD 221876 

5865 

4.29 

-0.60 

1.20 

-0.07 

-0.07 

0.05 

-0.03 + 0.07 

D 

L 






0.00 

-0.04 

0.11 

0.02 + 0.08 


N 

HD 224930 

5480 

4.45 

-0.66 

0.90 

-0.16 

-0.16 

-0.15 

-0.16 + 0.01 

7 

L 






-0.12 

-0.15 

-0.12 

-0.13 + 0.02 


N 

G69 - 8 

5640 

4.43 

-0.55 

1.10 

0.13 

0.07 

0.11 

0.10 + 0.03 

T 

L 






0.18 

0.09 

0.15 

0.14 + 0.05 


N 



